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Abstract 

Assume  a  constellation  of  satellites  is  flying  near  a  given  nominal  trajectory  around  L4  or  L5  in  the  Earth-Moon  system  in  such  a  way 
that  there  is  some  freedom  in  the  selection  of  the  geometry  of  the  constellation.  We  are  interested  in  avoiding  large  variations  of  the 
mutual  distances  between  spacecraft.  In  this  case,  the  existence  of  regions  of  zero  and  minimum  relative  radial  acceleration  with  respect 
to  the  nominal  trajectory  will  prevent  from  the  expansion  or  contraction  of  the  constellation.  In  the  other  case,  the  existence  of  regions  of 
maximum  relative  radial  acceleration  with  respect  to  the  nominal  trajectory  will  produce  a  larger  expansion  and  contraction  of  the  con¬ 
stellation.  The  goal  of  this  paper  is  to  study  these  regions  in  the  scenario  of  the  Circular  Restricted  Three  Body  Problem  by  means  of  a 
linearization  of  the  equations  of  motion  relative  to  the  periodic  orbits  around  L4  or  L5.  This  study  corresponds  to  a  preliminar  planar 
formation  flight  dynamics  about  triangular  libration  points  in  the  Earth-Moon  system.  Additionally,  the  cost  estimate  to  maintain  the 
constellation  in  the  regions  of  zero  and  minimum  relative  radial  acceleration  or  keeping  a  rigid  configuration  is  computed  with  the  use  of 
the  residual  acceleration  concept.  At  the  end,  the  results  are  compared  with  the  dynamical  behavior  of  the  deviation  of  the  constellation 
from  a  periodic  orbit. 
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1.  Introduction 

The  concept  of  Satellite  Formation  Flying  (SFF)  means 
to  have  two  or  more  satellites  in  orbit  such  that  their 
relative  positions  remain  constant  or  obeying  a  certain 
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dynamical  configuration  along  the  trajectory  (Sholomitslcy 
et  al.,  1977;  Battrick,  2000;  Bristow  et  ah,  2000;  Burns 
et  al.,  2000;  Ticker  and  Azzolini,  2000;  Fridlund  and 
Capaccioni,  2002).  This  concept  involves  the  control  over 
the  coordinated  motion  of  a  group  of  satellites,  with  the 
goal  of  maintaining  a  specific  geometric  space  configura¬ 
tion  between  the  elements  of  the  cluster  (Sabol  et  al., 
2001).  It  allows  that  a  group  of  low  cost  small  satellites, 
arranged  in  a  space  formation  flying,  operate  like  a  large 
Virtual  satellite’.  This  formation  will  have  many  benefits 
over  single  satellites,  including  simpler  designs,  faster  build 
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times,  cheaper  and  unprecedented  high  resolution  (Kapila 
et  ah,  2000).  Over  the  past  decade,  numerous  formation  fly¬ 
ing  missions  have  been  conceived.  The  Laser  Interferome¬ 
ter  Space  Antenna  (LISA)  is  a  proposed  mission,  whose 
objective  is  to  observe  astrophysical  and  cosmological 
sources  of  gravitational  waves  of  low  frequencies,  goal  that 
could  be  possible  using  three  identical  spacecraft  flying  in  a 
triangular  constellation,  with  equal  arms  of  5  million 
kilometers  each  (Peterseim  et  al.,  2000).  Another  example 
is  the  PRISMA  formation  flying  and  rendezvous  technol¬ 
ogy  mission,  launched  successfully  on  June  25,  2010,  into 
a  Sun  synchronous  orbit  at  approximately  750  km  altitude 
(Persson  et  al.,  2006;  Gill  et  al.,  2007;  Heilman  et  al.,  2009; 
Persson  et  al.,  2010).  PRISMA  mission  consists  of  two 
spacecraft:  Mango  and  Tango,  with  a  total  mass  of  about 
200  kg,  and  its  primary  purpose  is  to  demonstrate  forma¬ 
tion  flying  and  rendezvous  technology,  not  only  in  terms 
of  Guidance,  Navigation  and  Control  software  and  algo¬ 
rithms,  but  also  in  terms  of  instruments  and  operational 
aspects  (e.g.,  small  rocket  engines  and  Micro  Electro 
Mechanical  Systems).  Another  interesting  formation  flying 
mission  is  the  New  Worlds  Observer  (NWO)  (Cash  et  al., 
2009).  NWO  consists  of  a  large  telescope  and  an  occulter 
spacecraft  in  tandem  at  about  50,000  km  apart.  The  two 
spacecraft  would  be  flying  at  the  Earth-Sun  L2  Lagrangian 
point  or  in  a  drift-away  solar  orbit.  Its  purpose  is  to 
discover  and  analyze  terrestrial  extra-solar  planets.  The 
NWO  planned  launch  date  is  about  2018. 

A  configuration  of  SFF  can  typically  be  positioned  and 
maintained  in  two  dynamically  distinct  scenarios:  in  a  plan¬ 
etary  orbit  or  in  outer  space  (Alfriend  et  al.,  2002).  In  the 
planetary  orbit  scenario  the  fundamental  model  is  the 
problem  involving  two  light  bodies  (satellites),  close  to  each 
other,  that  could  be  in  the  same  orbit  or  describe  orbits  of 
slightly  different  radii,  eccentricities  and  inclinations 
around  a  massive  central  body,  e.g.  the  Earth.  Because 
the  gravitational  attraction  between  the  two  satellites  is 
negligible,  this  scenario  can  be  considered  like  a  superposi¬ 
tion  of  two  problems  of  two  bodies  (Sengupta  and  Vadali, 
2007),  where  the  atmospheric  drag  and  higher-order  grav¬ 
ity  terms  (e.g.,  J2)  can  be  considered  in  the  dynamic  models 
of  the  formation  (Carter  and  Humi,  2002;  Humi  and 
Carter,  2006;  Sengupta  et  al.,  2007)  and,  therefore,  improv¬ 
ing  the  fidelity  of  the  model.  An  example  of  this  kind  of 
formation  is  the  pair  of  satellites  Landsat  7  with  EO-1,  mis¬ 
sion  designed  to  enable  the  development  of  future  Earth 
imaging  observatories  that  will  have  a  significant  increase 
in  performance  while  reducing  cost  and  mass  (Flick,  2012). 

In  the  outer  space  scenario,  involving  SFF  in  the  vicinity 
of  a  libration  point,  the  fundamental  model  to  study  the 
natural  motion  of  spacecraft  relative  to  each  other  is  the 
Circular  Restricted  Three  Body  Problem  (CRTBP) 
(Hsiao  and  Scheeres,  2002).  The  interest  of  astronomical 
missions  has  been  to  position  SFF  around  the  Lagrangian 
points  L\  and  L2  (Faquhar,  1968;  Henon,  1973;  Breakwell 
and  Brown,  1979;  Howell,  1984;  Gomez  et  al.,  2000)  or 
L4  and  L5.  Of  particular  interest  is  the  stability  of  the  five 


Lagrangian  points.  In  the  case  of  the  three  collinear  sta¬ 
tionary  points:  L\,  L2  and  Z3,  they  are  always  unstable. 
Whereas  the  stability  of  L4  and  L5  points  depends  on  the 
mass  ratio  between  the  two  larger  bodies  (Danby,  1962; 
Szebehely,  1967).  In  important  cases  like  Earth-Moon  or 
Sun-Earth,  these  points  are  linearly  stable.  Moreover, 
there  exists  a  family  of  periodic  orbits  around  L4  and  L5. 
This  stability  property  makes  the  fuel  required  for  a  space¬ 
craft  to  maintain  its  relative  position  there  to  be  almost 
zero.  Despite  this  advantage,  today  there  are  no  missions 
orbiting  L4  or  L5  points  for  any  celestial  pair  of  primaries. 
In  the  case  of  Earth-Moon  system,  L4  and  L5  points  could 
be  excellent  locations  to  place  space  telescopes  for  astro¬ 
nomical  observations  or  a  space  station  (Schutz  et  al., 
1977).  In  addition,  there  is  a  renewed  interest  of  major 
space  agencies  for  Lagrangian  point  colonization.  Further¬ 
more,  Defilippi  (1977)  has  made  a  review  of  the  ideas  of 
O’Neill  (1974)  about  building  space  colonies  at  the  L4 
and  L5  positions.  These  space  stations  could  be  used  as  a 
way-point  for  traveling  to  and  from  the  region  between 
Earth’s  atmosphere  and  the  Moon  (cis-lunar  space). 

For  the  reason  that  keeping  a  formation  from  drifting 
apart  and  achieving  mission  requirements  is  expected  to 
require  significantly  more  fuel  than  station  keeping  a  single 
spacecraft,  one  of  the  problems  of  positioning  satellites  in 
formation  flying  is  the  cost  to  maintain  them  continuously 
orbiting  each  other.  In  particular,  missions  in  the  vicinity 
of  the  Lagrangian  points,  considering  the  scenario  of  the 
CRTBP,  may  be  placed  in  families  of  halo  orbits.  All  these 
orbits  are  inherently  unstable  (Gomez  et  al.,  1993)  and 
drive  the  SFF  out  of  its  desired  configuration.  Thereby,  a 
less  difficult  option  is  to  place  the  Satellite  Formation  Fly¬ 
ing  in  the  vicinity  of  Earth-Moon  L4  and  L5  Lagrangian 
points.  Because  of  their  stability  properties,  less  fuel  would 
be  spent  to  keep  the  formation  in  its  proper  configuration. 

Previous  studies  like  those  by  Catlin  and  McLaughlin 
(2007)  and  Wong  (2009)  on  SFF  about  L4  in  the  Earth- 
Moon  system  have  been  analyzed.  The  motion  of  formation 
flying  near  triangular  libration  points  was  studied  adopting 
the  CRTBP  model  and  the  linearized  equations  at  the 
equilibrium  points.  Catlin  and  McLaughlin  (2007)  show  that 
formations  are  possible  at  the  Triangular  points  on 
uncontrolled  trajectories  due  to  the  stability  of  stationary 
solutions.  On  the  other  hand,  Wong  (2009)  establishes  the 
need  for  a  system  control  and  develops  strategies  for  control¬ 
ling  a  spacecraft  formation  system  at  the  L4.  In  the  CRTBP 
scenario,  these  studies  show  that  velocity  change  require¬ 
ments  demanded  by  the  control  methods  would  be  very 
small.  Thus,  Catlin  and  McLaughlin  (2007)  and  Wong 
(2009)  conclude  that  nonlinear  aspects  as  well  as  perturba¬ 
tions  forces  (e.g.  solar  gravity,  solar  radiation  pressure)  are 
necessary  to  provide  a  more  real-world  accurate  descriptions 
of  motion  and  control  of  formation  flying  around  equilateral 
equilibrium  points  in  the  Earth-Moon  system. 

In  this  work,  we  do  not  adopt  the  linear  model  about  L4 , 
we  investigate  the  relative  dynamical  behavior  of  satellites 
with  respect  to  periodic  orbits  around  L4  and  L5  in  the 
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Earth-Moon  system  and  look  for  stable  regions  with  good 
properties  for  formation  flying.  In  particular,  we  are  inter¬ 
ested  in  finding  relative  positions  of  the  formation  with 
respect  to  the  periodic  orbits  such  that  large  variations  of 
the  mutual  distance  between  satellites  are  avoided.  In  this 
manner,  supposing  that  the  relative  velocity  of  each  satel¬ 
lite  with  respect  to  the  periodic  orbit  is  zero,  the  regions 
with  good  properties  for  formation  flying  are  such  that 
the  relative  radial  acceleration  with  respect  to  the  periodic 
orbit  is  zero  (Gomez  et  al.,  2006).  We  present  analytical 
and  numerical  methods  based  on  the  linearization  of  the 
relative  equations  of  motion  with  respect  to  periodic  orbits 
to  find  the  regions  where  the  radial  component  of  the  rela¬ 
tive  acceleration  is  zero.  Gomez  et  al.  (2006)  applied  a  sim¬ 
ilar  methodology  for  halo  orbits  around  the  collinear 
equilibrium  point  L2  in  the  Sun-Earth  system,  which  are 
unstable,  and  showing  the  existence  of  regions  of  Zero 
Relative  Radial  Acceleration  for  halo  orbits.  In  our  case, 
we  are  interested  in  investigating  the  existence  of  these 
regions  for  a  planar  formation  flight  dynamics  near  periodic 
orbits  around  the  equilateral  equilibrium  points  L4  and  L5  in 
the  Earth-Moon  system.  Additionally,  the  cost  to  maintain 
a  formation  in  these  regions  or  keeping  a  rigid  configuration 
is  analyzed  according  to  the  residual  acceleration  concept  as 
well  as  the  expansion  or  contraction  of  the  configuration.  In 
this  study,  we  consider  the  scenario  of  CRTBP  assuming 
that  all  spacecraft’s  trajectories  lie  on  the  same  orbital  plane, 
such  that,  Moon’s  eccentricity  and  perturbations  from  Sun 
will  be  neglected.  This  simplified  model  offers  a  good  first 
approximation  for  the  planar  formation  flight  dynamics 
about  triangular  libration  points.  However,  since  perturba¬ 
tions  from  Moon  and  Sun  affect  the  stability  of  triangular 
libration  points  (Schechter,  1968;  Kamel  and  Breakwell, 
1970;  Tapley  and  Schultz,  1970;  Gomez  et  al.,  1987; 
Gomez  et  al.,  2001;  McLaughlin  et  al.,  2004),  they  will  be 
considered  in  subsequent  works,  obtaining  a  more  accurate 
model,  which  is  important  for  practical  purposes. 

This  article  proceeds  as  follows:  In  Section  2,  we 
describe  the  equations  of  motion  of  the  formation  on  the 
plane  close  to  periodic  orbits  around  L4  and  L5 ;  next,  in 
Section  3,  we  describe  the  analytical  and  numerical 
methods  to  determine  the  regions  of  Zero  Relative  Radial 
Acceleration  as  well  as  we  define  mathematically  the 
residual  acceleration  concept.  In  Section  4,  we  determine 
the  regions  of  Zero,  minimum  and  maximum  relative  radial 
acceleration  for  two  kinds  of  periodic  orbit  families.  In 
Sections  5  and  6,  we  estimate  the  cost  to  maintain  a  space¬ 
craft  on  these  regions  and  keeping  a  rigid  configuration, 
respectively.  In  Section  7,  we  study  the  dynamical  behavior 
of  the  distance  between  the  spacecraft  and  the  periodic 
orbits  for  different  initial  conditions;  and  finally,  in 
Section  8,  we  present  our  conclusions  on  the  study. 

2.  Equations  of  local  dynamics  of  satellite  formation  flying 

In  this  section  we  describe  a  linear  approach  to  the  prob¬ 
lem  of  formation  flying  that  gives  the  relevant  information 


about  the  local  dynamics  of  the  problem.  In  order  to  avoid 
expansion  or  contraction  in  a  constellation  of  spacecraft, 
we  have  studied  the  existence  of  regions  with  Zero  Relative 
Radial  Acceleration  (ZRRA).  For  a  simple  model,  such  as 
the  CRTBP,  it  is  possible  to  compute  an  analytical  expres¬ 
sion  for  these  regions.  If  the  radius  of  the  constellation 
(largest  separation  between  spacecraft)  is  small,  then  a  lin¬ 
ear  approach  gives  all  the  relevant  information  about  the 
local  dynamics  of  the  problem  (Gomez  et  al.,  2006).  As  a 
first  approximation  to  study  the  natural  motion  of  space¬ 
craft  relative  to  each  other  near  periodic  orbits  around 
L4,  our  model  only  considers  motions  on  the  plane  of  the 
orbit  of  the  Earth-Moon  system.  As  reference  solutions, 
we  will  use  families  of  periodic  orbits  around  L4  since  the 
results  obtained  are  the  same  for  periodic  orbits  around 
L5  due  to  the  Mirror  Image  Theorem  applied  in  the 
Earth-Moon  system  (Miele,  1960).  Considering  the  sum 
of  the  masses  of  the  two  primaries,  the  distance  between 
them  and  also  making  the  Newton’s  gravitational  constant 
equal  the  one,  in  a  frame  rotating  with  the  primaries,  the 
equations  of  motion  for  a  spacecraft  in  the  non-dimen¬ 
sional  variables  are: 


..  ..  dU 
..  „ .  dU 

y  +  1'  =  w 

where  U  =  U(x,y)  is  given  by 

U(x,y)  =  -  {x2  +y2)  H - -  +  —  +  ~  p(l  -  it) 

2  Pi  p2  2 


(1) 

(2) 

(3) 


and  the  non-dimensional  relative  distance  pt  is  defined  as 

Pi  =  -Xif+y2.  (4) 

Let  A  =  (x,y,x,y)T  be  the  vector  that  describes  the  posi¬ 
tion  and  velocity  of  the  mass  m.  Hence  the  equations  of 
motion  (1)  and  (2)  can  be  written  in  vector  form  as 


X=f(X), 

where 

f  x  \ 

y 

f=  2y  +  f 

v-^+fy 


(5) 


(6) 


Now  assume  a  formation  such  that  the  chief  satellite  is 
placed  in  a  periodic  orbit  around  L4  as  show  in  Fig.  1. 
Let  Xh(t)  and X(t)  be  the  vectors  corresponding  to  the  tra¬ 
jectories  of  the  chief  and  deputy  satellites,  respectively. 
Then  the  vector  AX(t)  =  X(t)  —Xh(t)  represents  the  rela¬ 
tive  position  and  velocity  vectors  of  the  deputy  satellite 
with  respect  to  the  chief  satellite  at  instant  t.  The  coordi¬ 
nates  of  vector  AX(t)  are  defined  in  a  coordinate  system 
x',  /  parallel  to  the  rotating  coordinate  system  x,  y  and 
centered  at  Xh(t)  (see  Fig.  1). 
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Fig.  1.  Illustration  of  a  satellite  formation  flying  about  L4.  Xh(t)  and  X(t) 
denote  the  corresponding  trajectories  of  the  chief  and  deputy  satellites, 
respectively.  A X{t)  represents  the  relative  position  and  velocity  of  the 
deputy  satellite  at  instant  t  with  respect  to  a  coordinate  system  x,  y 
parallel  to  the  rotating  coordinate  system  x,  y  and  centered  at  Xh(t). 


Differentiating  the  vector  AX  with  respect  to  time  t ,  we 
have 

A X=X-Xh=f{X)-f(Xh).  (7) 

Supposing  that  the  diameter  of  the  formation  (largest 
separation  between  the  spacecraft)  is  small  compared  with 
the  chiefs  distances  with  respect  to  the  Earth  and  the 
Moon,  px  and  p2 ,  respectively,  then  we  can  linearize  the 
vector  field  /(X)  around  the  periodic  orbit  Xh{t)  (Hill, 
1878;  Clohessy  and  Wiltshire,  1960) 

f(X)nf(Xh)+Df(Xh)  AX.  (8) 

Substituting  this  expression  in  Eq.  (7),  we  obtain  the  lin¬ 
ear  behavior  of  the  deputy  satellite  around  a  periodic 
solution: 


AX  =  Df(Xh)AX, 

where 

(  ° 

0 

l 

°\ 

0 

0 

0 

1 

Df  = 

d2U 

dx2 

d2u 

dydx 

0 

2 

d2U 
\  dxdy 

d2U 

dy2 

-2 

0  / 

(9) 


(10) 


3.  Zero  Relative  Radial  Accelerations  Lines 


d2U 

d2U 

dx2 

dydx 

d2U 

d2u 

dxdy 

dy2 

The  points  with  zero  relative  velocity  are  those  such  that 
v  =  0,  and,  in  this  case,  we  have  that  the  relative  accelera¬ 
tion  is  given  by 

r  =  Fr.  (12) 

Therefore,  the  radial  component  of  the  relative  acceler¬ 
ation  will  be  zero  in  the  set  of  points  where  the  vectors  f 
and  r  are  orthogonal,  in  other  words, 

rTFr  =  0.  (13) 


Eq.  (13)  represents  two  lines  which  depend  on  the  point 
Xh(t )  selected  along  the  periodic  solution  of  system  (5). 

Zero  Relative  Radial  Acceleration  Lines  (ZRRAL)  can 
also  be  computed  numerically  (Gomez  et  al.,  2006).  Given 
a  certain  periodic  orbit,  we  select  a  point  on  it:  Xh(t)  = 
{x/i{t),vh(t))T .  Around  this  point,  we  consider  a  sphere  in 
the  configuration  space,  of  radius  s  and  we  set  the  velocity 
of  all  the  points  of  the  sphere  equal  to  the  velocity  of  the 
chosen  point,  Vh(t)  (zero  relative  velocity  condition).  Using 
polar  coordinates,  the  set  of  points  on  the  sphere  are  given 
by 


/  Xh(t )  +  s(0)  \ 

v  vk(t)  y 


(14) 


Now,  writing  the  equations  of  the  motion  (5)  for  Xh  as 


(xh\  _  ffi{xh,vh)\ 
\xh)  \f2{xh,vh)J' 


(15) 


the  relative  acceleration  can  be  evaluated  by  means  of 
*(t,  0)  =f2(xh(t),  vh{t))  -f2{xh{t)  +  s(0),vh(t)),  (16) 


whose  scalar  product  with  s(0)  will  be  the  desired  relative 
radial  acceleration  for  each  angle  6.  In  this  manner,  we 
denote  by  6 *  the  angle  at  which  the  scalar  product  between 
a(t,  6 *)  and  s(0*)  is  equal  to  zero  and,  therefore,  it  will  give 
us  the  relative  position  vector  v  that  belongs  to  ZRRAL. 

Now,  we  will  show  how  to  estimate  the  cost  to  maintain 
spacecraft  in  formation  using  the  concept  of  residual 
acceleration.  This  concept  can  be  used  not  only  if  the 
formation  follows  the  direction  of  ZRRAL  but  for  any 
rigid  configuration. 


3.1.  Residual  acceleration  and  cost  estimate  to  maintain  a 
spacecraft  in  a  formation 


Assume  two  satellites  flying  in  a  formation  in  a  periodic 
orbit  around  Z4.  Writing  the  array  AX  as  (r,  v)T,  where 
r  =  (Ax,  Ay)r,  the  linear  system  (9)  becomes 


) 


0  / 
F  J 


(11) 


where 


Given  a  certain  nominal  trajectory  Xh(t)  around  Z4,  sup¬ 
pose  that  a  spacecraft  follows  an  artificial  trajectory  around 
it,  i.e.  there  exists  a  continuous  effort  (control)  applied  on  it 
necessary  to  maintain  that  trajectory.  Denote  by  ra(t)  the 
relative  position  at  time  t.  If  there  is  no  continuous  effort 
applied  on  the  spacecraft,  then  the  function  ra(t)  will  satisfy 
Eq.  (11).  However,  if  a  control  is  applied  on  the  spacecraft, 
there  is  a  residual  acceleration  as  a  consequence  of  the 
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maneuvers  applied  on  it.  In  this  manner,  denoting  by  Ra  the 
residual  acceleration,  it  can  be  computed  using  Eq.  (11): 

Ra  =  ra-  [Fra  +  Jra).  (17) 


Therefore,  the  cost  estimate  to  maintain  the  spacecraft 
on  an  artificial  trajectory  until  a  time  T  is  simply: 


AV  = 


(18) 


where  ||  ||  denotes  the  Euclidean  norm. 

In  this  simple  manner,  we  can  evaluate  the  cost  estimate 
to  maintain  two  or  more  satellites  in  formation  along  a 
periodic  orbit  around  Z4. 


4.  Zero  Relative  Radial  Acceleration  Lines  in  a  long  and 
short  period  family 

In  this  section  we  have  determined  the  existence  of  Zero 
Relative  Radial  Acceleration  Lines  (ZRRAL)  in  periodic 
orbits  around  L4.  The  existence  of  ZRRAL  along  any  nom¬ 
inal  trajectory  is  determined  by  the  sign  of  the  discriminant 
of  the  sub-matrix  F  in  Eq.  (13)  which  represents,  in  general, 
a  quadratic  form.  The  discriminant  of  F  is  defined  by  the 
following  expression:  B2  —  4 AC,  where  A  =  d2U/dx 2,  B  = 
ld2U /dxdy,  and  C  =  d2U/dy2.  If  the  discriminant  of  F  is 
negative  at  a  certain  point  of  the  nominal  trajectory,  then, 
there  is  no  region  with  Zero  Relative  Radial  Acceleration 
(ZRRA).  Otherwise,  the  ZRRA  at  this  point  is  represented 
by  two  lines  (in  the  planar  case)  (Lawrence,  1972;  Barry, 
2007). 

First,  we  begin  by  computing  the  value  of  the  discrimi¬ 
nant  of  the  sub-matrix  F  associated  to  points  Xh{t)  in 
two  periodic  orbits  close  enough  to  Z4,  as  shown  in 
Fig.  2.  Their  period  is  92  days  (long  period,  Fig.  2(a)) 
and  28  days  (short  period,  Fig.  2(b)).  These  two  orbits,  as 
well  as  all  of  the  periodic  orbits  shown  in  this  section,  are 
seen  with  respect  to  a  non-dimensional  coordinate  system 


(a) 


(b) 


Fig.  2.  Elliptic  orbits  centered  at  L4  of  (a)  92  days  and  (b)  28  days  period. 
They  correspond  to  the  well  known  families  of  long  and  short  period 
orbits  about  L4  in  the  CRTBP. 


x",  y"  rotated  29.7°  clockwise  about  the  equilibrium  point 
U  (Szebehely,  1967). 

As  we  can  see  in  Fig.  3,  the  sign  of  the  discriminant  in 
both  trajectories  are  negative  for  every  point  Xh(t), 
meaning  that  there  are  no  regions  with  ZRRA  along  these 
specific  periodic  orbits. 

The  previous  fact  can  be  verified  when  we  compute 
numerically  the  ZRRA  at  some  points  on  these  periodic 
orbits.  Thus  considering  a  sphere  of  radius  s  equal  to 
1  km,  the  first  row  of  Fig.  4  shows  the  scalar  product 
between  the  relative  acceleration  a(t,  6)  and  vector  s(0)  as 
a  function  of  angle  6  at  three  different  points  Xh(t)  along 
the  long  period  orbit,  where  t  =  22,  45  and  68  days. 
Similarly,  the  second  row  of  Fig.  4  shows  the  scalar 
product  at  three  different  points  Xh(t)  along  the  short 
period  orbit,  where  t  =  7,  14  and  22  days. 

Now,  if  we  make  a  zoom  around  the  maximum  point  of 
each  dot  product  function  shown  in  Fig.  4,  the  scalar  prod¬ 
uct  function  between  the  relative  acceleration  and  relative 
position  vectors  never  crosses  the  horizontal  axis  for  any 
angle  6  at  each  of  these  three  points  Xh(t),  as  we  can  see 
in  Fig.  5.  Therefore  the  radial  component  of  the  relative 
acceleration  a(t,  6)  is  different  from  zero  for  any  angle  6 
at  each  one  of  these  three  points.  It  is  important  to  com¬ 
ment  that  the  dot  product  computed  in  this  work  is  dimen¬ 
sionless  because  we  are  using  the  normalized  system.  If  we 
want  to  obtain  the  scalar  product  in  physical  units,  i.e. 
meters  and  seconds  (time  and  length),  we  have  to  multiply 
the  previous  results  by  a  factor  of  1.0469  x  106  m/s2.  In 
any  case,  the  qualitative  behavior  of  the  scalar  product 
function  shown  in  Fig.  4  will  remain  the  same.  Addition¬ 
ally,  from  Fig.  4  we  can  see  that  the  scalar  product  function 
is  periodic  (period  equal  to  180°)  and  has  two  maxima  and 
two  minima.  Thus,  the  radial  component  of  the  relative 
acceleration  a(t,  6)  has  also  two  maxima  and  two  minima 
at  the  same  previous  optimal  values  of  0.  In  Fig.  6  we  show 
the  angle  0  for  each  point  Xh  ( t )  where  the  relative  radial 
acceleration  has  a  maximum  and  a  minimum  value  along 
the  long  (first  row)  and  short  (second  row)  period  orbits 
shown  in  Fig.  2. 

From  this  first  case,  we  could  state  that  there  are  no 
regions  with  ZRRA  along  the  periodic  orbits  that  are  close 
enough  to  Z4.  The  next  step,  therefore,  is  to  explore  the 
existence  of  regions  with  ZRRA  along  periodic  orbits  that 
are  far  from  Z4. 

Fig.  7  shows  two  periodic  orbits  around  Z4,  which 
belong  to  the  long  (92  days)  and  short  (21  days)  period 
families,  in  the  non-dimensional  coordinate  system  x",  y" . 
Both  are  much  larger,  i.e.  the  amplitude  is  larger  when 
compared  to  the  periodic  orbits  shown  in  Fig.  2. 

Similarly,  Fig.  8(a)  and  (b)  show  the  value  of  the  dis¬ 
criminant  of  the  sub-matrix  F  associated  to  the  points 
Xh(t)  along  the  long  and  short  period  orbits  shown  in 
Fig.  7.  Unlike  the  previous  case,  the  sign  of  the  discrimi¬ 
nant  in  both  trajectories  are  negative  and  positive,  this 
means  that  although  there  do  not  exist  regions  with  ZRRA 
at  certain  points  along  the  periodic  orbits,  there  exist  some 
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Fig.  3.  Discriminant  of  the  sub-matrix  F  associated  to  the  points  Xh(t)  in  the  (a)  long  and  (b)  short  period  orbits  shown  in  Fig.  2. 


0  (deg)  0  (deg)  0  (deg) 

Fig.  4.  Scalar  product  between  the  relative  acceleration  a(t,  0)  and  vector  s(6)  as  a  function  of  angle  0  at  three  different  points  Xh(t)  along  the  long  (first 
row)  and  short  (second  row)  period  orbits  shown  in  Fig.  2. 


-i3  t  =  22  days 


-i3  t  =  68  days 


0-i3  t  =  22  days 


0  (deg) 


Fig.  5.  Scalar  product  between  the  relative  acceleration  a(t,  0)  and  vector  s(6)  as  a  function  of  angle  0  at  three  different  points  Xh(t)  along  the  long  (first 
row)  and  short  (second  row)  period  orbits  shown  in  Fig.  2. 
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Maximum  Minimum 


Fig.  6.  Angle  0  for  each  point  Xh(t)  where  the  relative  radial  acceleration  has  a  maximum  and  minimum  value  on  the  long  (first  row)  and  short  (second 
row)  period  orbits  shown  in  Fig.  2. 


parts  along  these  periodic  orbits  where  the  radial  compo¬ 
nent  of  the  relative  acceleration  is  equal  to  zero.  This  fact 
can  also  be  seen  if  we  consider  a  sphere  of  radius  s  equal 
to  1  km  as  shown  in  Fig.  9  and  we  compute  numerically 


the  ZRRA  at  three  different  points Xh(t)  along  the  previous 
periodic  orbits. 

The  qualitative  behavior  of  the  scalar  product  function 
is  very  similar  for  all  the  values  of  t  where  there  appear 
two  maxima  and  two  minima.  However,  there  is  a  set  of 
points  in  these  trajectories  where  the  radial  component  of 
the  function  a{t ,  6)  is  zero  with  vertex  at  Xh(t )  for  two  dif¬ 
ferent  values  of  0 ,  which  we  have  denoted  by  6 *  and  6 ** 
(where  0*  ^  0**).  Note  that  since  the  scalar  product  func¬ 
tion  is  periodic  with  respect  to  6 ,  the  other  two  zeros  rep¬ 
resent  the  same  situation.  Fig.  10(a)  and  (b)  show  the 
function  0  associated  to  point  Xh(t),  such  that,  the  relative 
radial  acceleration  is  maximum  for  all  values  of  t  along  the 
long  and  short  period  orbits  shown  in  Fig.  7,  respectively. 
Similarly,  Fig.  1 1  shows  the  function  6  for  the  values  where 
the  relative  radial  acceleration  is  minimum  or  zero  for  the 
long  (first  row)  and  short  (second  row)  period  orbits.  As  we 
can  see  in  Fig.  11,  when  we  only  compute  the  values  of  0 
where  the  scalar  product  function  has  a  minimum,  the 
function  6  is  smooth.  But,  when  we  compute  the  points 
where  the  scalar  product  function  has  a  minimum  or  a 
zero,  there  are  two  points  where  this  function  is  not 


O  lO  20 

t  (days) 

(b) 


Fig.  8.  Discriminant  of  the  sub-matrix  F  associated  to  the  points  in  the  (a)  long  and  (b)  short  period  orbits  shown  in  Fig.  7. 
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Fig.  9.  Scalar  product  between  the  relative  acceleration  a (t,  0)  and  vector  s (0)  as  a  function  of  angle  6  at  three  different  points Xh(t)  on  the  long  (first  row) 
and  short  (second  row)  period  orbits  shown  in  Fig.  7. 


t  (ttlivs) 
(b) 


Fig.  10.  Angle  0  for  each  point  Xh(t)  where  the  relative  radial  acceleration  has  a  maximum  value  for  the  (a)  long  and  (b)  short  period  orbits  shown  in 
Fig.  7. 


Minimum  Minimum  and  0*  Minimum  and  0** 


Minimum  Minimum  and  0*  Minimum  and  0** 


Fig.  11.  Angle  0  for  each  point  Xh(t)  where  the  relative  radial  acceleration  has  a  minimum  or  zero  value  for  the  long  (first  row)  and  short  (second  row) 
period  orbits  shown  in  Fig.  7. 
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smooth.  This  fact  implies  a  higher  cost  to  maintain  a  rigid 
constellation  as  will  be  shown  later. 

Finally,  Fig.  12(a)  and  (b)  show  the  family  of  long  and 
short  period  around  Z4.  Similarly,  we  compute  the  discrim¬ 
inant  of  the  sub-matrix  F  along  each  periodic  orbit  as 
shown  in  Fig.  13.  This  result  supports  a  fact  that  was 


(b) 

Fig.  12.  (a)  Long  and  (b)  short  period  families  around  Z4.  Their  periods 
are  about  92  days  and  28  days,  respectively. 


previously  shown.  When  their  amplitude  is  small  enough, 
there  are  no  ZRRAL.  However,  when  their  amplitude  is 
large  enough,  ZRRA  regions  exist  and  their  size  increases 
with  the  amplitude. 

We  have  computed  the  value  of  0  for  every  point  Xh(t) 
along  each  orbit  that  belongs  to  the  long  and  short  period 
families  where  the  radial  component  acceleration  is  maxi¬ 
mum,  minimum  or  zero  (in  the  cases  where  the  discrimi¬ 
nant  is  positive).  Fig.  14(a)  and  (b)  show  the  value  of  0 
such  that  the  radial  acceleration  component  is  maximum 
for  the  long  and  short  period  families,  respectively. 

In  the  same  way,  for  periodic  orbits  that  are  large 
enough,  the  function  9  has  a  minimum  and  two  zeros. 
Fig.  15  shows  the  value  of  6  such  that  the  radial  accelera¬ 
tion  component  is  minimum  or  zero  (denoted  by  6 *  and 
6 **)  for  the  long  (first  row)  and  short  (second  row)  period 
families,  respectively. 

5.  Cost  estimate  to  maintain  a  spacecraft  on  the  regions  of 
zero,  minimum  and  maximum  relative  radial  accelerations 

Considering  a  radius  s  equal  to  1  km,  Table  1  shows  the 
cost  estimate  per  year  to  maintain  the  spacecraft  when  the 
artificial  trajectory  follows  the  direction  where  the  relative 
radial  acceleration  component  is  maximum  (see  left  column 


Fig.  13.  Discriminant  of  the  sub-matrix  F  associated  to  the  points  in  the  (a)  long  and  (b)  short  period  families  shown  in  Fig.  12. 


t  (days) 
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Fig.  14.  Angle  6  for  each  point Xh(t)  where  the  relative  radial  acceleration  has  a  maximum  value  for  the  (a)  long  and  (b)  short  period  families  shown  in 
Fig.  12. 
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Fig.  15.  Angle  0  for  each  point  Xh(t)  where  the  relative  radial  acceleration  has  a  minimum  or  zero  value  for  the  long  (first  row)  and  short  (second  row) 
period  families  shown  in  Fig.  12. 


Table  1 

Cost  estimate  per  year  to  maintain  a  spacecraft  considering  a  radius  of 
1  km  and  following  the  direction  where  the  relative  radial  acceleration 
component  is  maximum  and  minimum  with  respect  to  the  periodic  orbits 
shown  in  Fig.  2. 


Orbital  period 

Cost  estimate 

Cost  estimate 

Maximum  radial 

Minimum  radial 

acceleration  (m/s) 

acceleration  (m/s) 

Long 

6.5777  x  10“‘ 

6.0422  x  10“3 

Short 

6.2989  x  10“‘ 

6.1713  x  10'3 

in  Fig.  6)  and  minimum  (see  right  column  in  Fig.  6)  with 
respect  to  the  periodic  orbits  shown  in  Fig.  2.  It  was  found 
that  there  exists  a  linear  relation  between  the  distance  from 
Xh(t )  and  the  cost  estimate  when  the  radius  s  is  small 
enough.  Additionally,  the  costs  between  the  long  and  short 
period  orbits  are  almost  equal. 

Now,  suppose  a  spacecraft  placed  in  an  artificial  trajec¬ 
tory  around  the  periodic  orbits  shown  in  Fig.  7.  Similarly, 
considering  a  radius  s  equal  to  1  km,  Table  2  shows  the  cost 
estimate  per  year  to  maintain  the  spacecraft  following  the 
direction  where  the  relative  radial  acceleration  component 
is  maximum  (see  Fig.  10),  minimum  and  zero  (see  Fig.  11). 
Similarly,  there  exists  a  linear  relation  between  the  cost  and 
the  radius  of  the  artificial  trajectory  when  y  is  small  enough. 


Note  that  the  cost  estimate  to  maintain  the  spacecraft,  on 
the  direction  when  the  relative  radial  acceleration  is  mini¬ 
mum  but  not  zero,  is  much  less  (almost  ten  times)  when 
compared  to  the  cost  to  maintain  the  spacecraft  on  the 
direction  where  the  relative  radial  acceleration  is  minimum 
or  zero  (see  Table  2).  This  fact  can  be  understood  if  we  com¬ 
pute  the  components  x  and  y  of  the  relative  position,  veloc¬ 
ity  and  acceleration  vectors  described  in  Eq.  (17).  For 
example,  Fig.  16  shows  the  components  x  (left  column) 
and  y  (right  column)  of  the  relative  position  va  (first  row), 
relative  velocity  ra  (second  row)  and  relative  acceleration 
va  (third  row)  of  the  artificial  trajectory  following  the  direc¬ 
tion  where  the  relative  radial  acceleration  component  is 
minimum  on  the  long  period  orbit  shown  in  Fig.  7.  Simi¬ 
larly,  Fig.  17  shows  the  components  of  these  vectors  when 
the  relative  radial  acceleration  component  is  minimum  or 
zero  on  the  long  period  orbit  shown  in  Fig.  7.  Note  that 
in  Fig.  16  the  two  components  for  rfl,  iflJ  ra ,  are  smooth 
functions.  However,  in  Fig.  17  we  can  see  that  the  compo¬ 
nents  for  each  vector  are  not.  In  fact,  the  function  ra  is 
not  smooth  at  time  t  when  the  regions  of  ZRRA  appear 
and  disappear  (see  Fig.  1 1).  Therefore,  we  see  two  big  jumps 
in  the  functions  va  and  ra  at  these  two  points.  The  function 
Ra  is  bigger  in  the  region  where  the  ZRRA  exists  as  shown 
in  Fig.  18.  As  a  consequence,  the  cost  estimate  AV  will  be 
less  than  that  when  the  artificial  trajectory  follows  the 


Table  2 


Cost  estimate  per  year  to  maintain  a  spacecraft  considering  a  radius  of  1  km  following  the  direction  where  the  relative  radial  acceleration  component  is 
maximum,  minimum  and  zero  with  respect  to  the  periodic  orbits  shown  in  Fig.  7. 


Orbital  period 

Cost  estimate 

Cost  estimate 

Cost  estimate 

Cost  estimate 

Maximum  radial 

Minimum  radial 

Minimum  and  zero  (0*) 

Minimum  and  zero  (0**) 

acceleration  (m/s) 

acceleration  (m/s) 

radial  acceleration  (m/s) 

radial  acceleration  (m/s) 

Long 

6.6074  x  10“‘ 

7.0033  x  10“3 

6.4182  x  10“2 

6.5106  x  10“2 

Short 

6.3015  x  IQ-1 

2.1255  x  10~2 

1.3563  x  10“' 

1.3635  x  IQ-1 
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Fig.  16.  Components  x  (right  column)  and  y  (left  column)  of  relative  position  ra  (first  row),  relative  velocity  ra  (second  row)  and  relative  acceleration  ra 
(third  row)  of  the  artificial  trajectory  that  follows  the  direction  where  the  relative  radial  acceleration  component  is  minimum  in  the  long  period  orbit  shown 
in  Fig.  7. 
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Fig.  17.  Components  x  (right  column)  and  y  (left  column)  of  relative  position  ra  (first  row),  relative  velocity  ra  (second  row)  and  relative  acceleration  ra 
(third  row)  of  the  artificial  trajectory  that  follows  the  direction  where  the  radial  acceleration  component  is  minimum  or  zero  in  the  long  period  orbit  shown 
in  Fig.  7. 


direction  where  the  relative  radial  acceleration  is  minimum 
but  not  zero. 

Similarly,  suppose  that  we  denote  the  family  of  periodic 
orbits  as  shown  in  Fig.  19,  where  Orbit  1  denotes  the  small¬ 
est  one  and  Orbit  8  denotes  the  largest  one.  Considering  a 
radius  equal  to  1  km,  Tables  3  and  4  show  the  cost  estimate 
per  year  to  maintain  a  spacecraft  in  each  orbit  of  the  long 
and  short  period  family,  respectively,  following  the 
direction  where  the  relative  radial  acceleration  component 
is  maximum,  minimum  and  zero.  Additionally,  the 


magnitude  of  the  segment  where  each  orbit  intersects  the 
positive  v-axis  is  also  shown.  This  quantity  is  useful  to  rep¬ 
resent  the  amplitude  of  each  orbit.  As  can  be  seen  in  Tables 
3  and  4,  the  cost  estimate  following  the  direction  where  the 
radial  acceleration  component  is  maximum  practically 
remains  constant,  i.e.,  it  does  not  depend  on  the  amplitude 
of  the  periodic  orbits.  On  the  other  hand,  following  the 
direction  where  the  radial  acceleration  component  is  mini¬ 
mum  or  zero,  the  cost  to  maintain  a  spacecraft  increases 
with  the  amplitude. 
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Fig.  18.  Residual  acceleration  of  the  artificial  trajectory  that  follows  the  direction  where  the  relative  radial  acceleration  component  is  minimum  or  zero 
with  respect  to  the  long  period  orbit  shown  in  Fig.  7. 


x " 

(b) 


Fig.  19.  Family  of  periodic  orbits  shown  in  Fig.  12,  where  Orbit  1  denotes 
the  smallest  one  and  Orbit  8  denotes  the  largest  one. 

Finally,  it  is  important  to  comment  that  although  all 
simulations  made  in  this  section  show  that  the  cost  estimate 
to  maintain  a  satellite  in  the  region  of  maximum  relative 
radial  acceleration  is  almost  ten  times  higher  than  the  cost 
estimate  to  maintain  it  in  the  region  of  minimum  relative 
radial  acceleration,  the  values  estimated  are  very  low  (the 
residual  acceleration  is  order  of  1(T8  m/s2),  such  small 


continuous  acceleration  to  control  the  relative  position  of 
a  satellite  can  not  be  implemented  in  practical  engineering 
applications.  In  this  manner,  effects  of  the  Moon’s  eccen¬ 
tricity  and  the  Sun  must  be  considered  in  the  dynamical 
model  to  correct  the  cost  estimate  to  maintain  the  forma¬ 
tion.  However,  as  will  be  shown  in  Section  7,  the  signifi¬ 
cance  of  the  ZRRAL  is  not  only  in  minimizing  the  cost 
of  correction  maneuvers  but  also  in  determining  the  initial 
conditions  of  the  SFF  that  minimize  the  drift  of  the  forma¬ 
tion.  In  practice,  it  is  difficult  to  find  initial  conditions  that 
eliminate  inter-spacecraft  drift  because  of  perturbations. 
By  contrast,  the  initial  conditions  can  be  picked  up  in  such 
a  way  minimizing  the  need  for  corrective  action. 

6.  Cost  estimate  to  maintain  a  formation  on  a  rigid 
configuration 

In  the  previous  section,  we  studied  the  cost  estimate  of 
configurations  that  follows  the  direction  where  the  relative 
radial  acceleration  is  maximum,  minimum  or  zero.  How¬ 
ever,  many  space  missions  (e.g.,  interferometry  missions) 
require  fixed  relative  positions  and  specific  geometric  tem¬ 
plates.  Therefore,  in  this  section  we  study  the  case  when  the 
geometry  of  the  configuration  remains  constant,  i.e.  the 
angle  6  is  fixed  (see  Fig.  24). 


Table  3 


Cost  estimate  per  year  to  maintain  a  spacecraft  in  each  orbit  of  the  long  period  family  shown  in  Fig.  19  considering  a  radius  equal  to  1  km  and  following 
the  direction  where  the  relative  radial  acceleration  component  is  maximum,  minimum  and  zero. 


Periodic  orbit 

Amplitude  (km) 

Cost  estimate 
Maximum  radial 
acceleration  (m/s) 

Cost  estimate 
Minimum  radial 
acceleration  (m/s) 

Cost  estimate 

Minimum  and  zero  (6*) 
radial  acceleration  (m/s) 

Cost  estimate 

Minimum  and  zero  (0**) 
radial  acceleration  (m/s) 

1 

3.6614  x  103 

6.5750  x  10“' 

6.0286  x  10“3 

2 

2.1526  x  104 

6.5792  x  10~‘ 

6.0881  x  10~3 

1.5004  x  10'2 

1.5062  x  10~2 

3 

3.8517  x  104 

6.5850  x  10  1 

6.2313  x  10~3 

3.3591  x  10~2 

3.3850  x  10~2 

4 

5.4585  x  104 

6.5925  x  10~‘ 

6.4562  x  10'3 

4.5180  x  10'2 

4.5567  x  10'2 

5 

6.9807  x  104 

6.6014  x  10“' 

6.7676  x  10“3 

5.5834  x  10“2 

5.6401  x  10“2 

6 

8.4300  x  104 

6.6116  x  10~‘ 

7.1855  x  10~3 

7.3302  x  10~2 

7.4046  x  10~2 

7 

9.8100  x  104 

6.6230  x  10“' 

7.7586  x  10“3 

8.3866  x  10“2 

8.4787  x  10“2 

8 

1.1171  x  105 

6.6360  x  10~‘ 

8.6312  x  10'3 

8.7558  x  10~2 

8.8503  x  10'2 
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Table  4 

Cost  estimate  per  year  to  maintain  a  spacecraft  in  each  orbit  of  the  short  period  family  shown  in  Fig.  19  considering  a  radius  equal  to  1  km  and  following 
the  direction  where  the  relative  radial  acceleration  component  is  maximum,  minimum  and  zero. 


Periodic  orbit 

Amplitude  (km) 

Cost  estimate 
Maximum  radial 
acceleration  (m/s) 

Cost  estimate 
Minimum  radial 
acceleration  (m/s) 

Cost  estimate 

Minimum  and  zero  ( 0 *) 
radial  acceleration  (m/s) 

Cost  estimate 

Minimum  and  zero  ( 0 **) 
radial  acceleration  (m/s) 

1 

2.1238  x  103 

6.4023  x  10'1 

5.9153  x  10~3 

2 

1.2570  x  104 

6.4013  x  10’1 

7.4009  x  10~3 

6.4286  x  10'2 

6.4920  x  10~2 

3 

2.2641  x  104 

6.4025  x  10“‘ 

1.0242  x  10“2 

9.1027  x  10“2 

9.1375  x  10“2 

4 

3.2443  x  104 

6.4054  x  10’1 

1.3850  x  10~2 

9.9301  x  10'2 

1.0013  x  10’1 

5 

4.1900  x  104 

6.4092  x  10“‘ 

1.7736  x  10“2 

1.2289  x  10_1 

1.2372  x  10_1 

6 

5.1048  x  104 

6.4137  x  10’1 

2.1740  x  10~2 

1.3883  x  10'1 

1.3987  x  10’1 

7 

5.9851  x  104 

6.4189  x  10“‘ 

2.5826  x  10“2 

1.4365  x  10“‘ 

1.4685  x  10_1 

8 

6.8347  x  104 

6.4255  x  10“’ 

2.9996  x  10~2 

1.8064  x  10'1 

1.8140  x  10*1 

Considering  a  radius  of  1  km,  Fig.  20(a)  and  (b)  show 
the  cost  estimate  per  year  to  maintain  a  spacecraft  in  each 
orbit  of  the  long  and  short  period  family  shown  in  Fig.  19, 
keeping  the  angle  0  constant. 

Note  that  the  cost  function  in  Fig.  20  has  a  period  of 
180°.  Additionally,  the  function  has  two  maxima  at 
6  «  60°  and  240°,  and  two  minima  at  6  «  150°  and  330°. 
These  values  practically  coincide  with  the  direction  in 
which  the  relative  radial  acceleration  component  is  maxi¬ 
mum  and  minimum  (see  Fig.  6)  for  the  periodic  orbits 
shown  in  Fig.  2.  In  fact,  the  costs  compared  with  these 
cases  are  similar. 

Now,  we  will  study  a  triangular  geometry  around  the 
periodic  orbits  shown  in  Figs.  2  and  7.  In  this  configuration 
we  consider  a  formation  of  three  satellites  such  that  their 
positions  and  distances  with  respect  to  periodic  orbit  are 
fixed  and  form  an  equilateral  triangle.  There  are  three  cases 
we  have  analyzed  in  this  formation  as  shown  in  Fig.  21: 


(a)  The  periodic  solution  Xh(t)  remains  in  the  center  of 
the  equilateral  triangle. 

(b)  The  periodic  solution  Xh{t)  remains  in  the  middle  of 
one  of  the  sides. 

(c)  One  of  the  satellites  remains  on  the  periodic  solution 

Xh(t). 

In  these  three  cases,  the  orientation  of  the  formation  is 
defined  by  the  angle  6 ,  which  denotes  the  direction  of  the 
relative  position  of  satellite  1.  Once  the  angle  0  is  fixed, 
the  other  angles  that  define  the  relative  positions  of  satel¬ 
lites  2  and  3  are  easily  computed  by  elementary  geometry. 
The  cost  estimate  to  maintain  the  formation  is  simply  the 
sum  of  the  cost  estimates  to  maintain  each  satellite  in  the 
selected  position.  Fig.  22(a)  and  (b)  show  the  cost  function 
of  the  triangular  formation  along  the  long  and  short  period 
orbits  around  L4  shown  in  Fig.  2,  for  0°  <  6  <  180°  (the 
cost  function  has  a  period  of  180°)  and  a  time  T  =  1  year. 


140  150  160  170 


0  (deg) 


Fig.  20.  Cost  estimate  per  year  to  maintain  a  spacecraft  in  the  long  (a)  and  short  (b)  period  family  shown  in  Fig.  12  fixing  the  angle  6  and  considering  a 
radius  of  1  km. 
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Fig.  21.  Formation  of  three  satellites  such  that  their  positions  and  distances  in  the  periodic  orbit  are  fixed  and  form  an  equilateral  triangle. 


(a)  (b) 


Fig.  22.  Cost  estimate  per  year  to  maintain  a  formation  of  three  satellites  considering  configurations  (a),  (b)  and  (c)  in  the  long  and  short  period  orbits 
shown  in  Fig.  2. 


The  length  of  the  side  of  the  equilateral  triangle  is  equal  to 
2  km.  As  we  can  see  in  Fig.  22,  for  the  long  and  short  per¬ 
iod  orbits,  the  minimum  cost  is  obtained  when  the  forma¬ 
tion  follows  the  configuration  (c)  for  6  =  90°  and  6  =  150°. 
However,  we  note  that  the  cost  of  configuration  (b)  is 
almost  the  same  when  6  =  6 0°.  Therefore,  it  can  be  said 


that  the  minimum  cost  is  obtained  when  the  formation  fol¬ 
lows  the  configuration  (b)  and  (c)  and  one  of  the  satellites  is 
located  along  the  direction  of  0  =  150°. 

Finally,  Fig.  23(a)  and  (b)  show  the  cost  estimate  func¬ 
tion  of  the  triangular  formation  on  the  long  and  short  per¬ 
iod  orbits  around  Z4  shown  in  Fig.  7.  The  results  show  that 


(a) 


144  145 

0  (deg) 

(b) 


Fig.  23.  Cost  estimate  per  year  to  maintain  a  formation  of  three  satellites  considering  configurations  (a),  (b)  and  (c)  in  the  long  and  short  period  orbits 
shown  in  Fig.  7. 
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the  minimum  cost  estimate  is  obtained  when  the  formation 
follows  the  configuration  (c)  for  6  =  103°  and  6  =  143°  for 
the  long  and  short  period  orbits,  respectively.  Additionally, 
it  is  interesting  to  remark  that  the  configuration  (c)  is  not 
only  the  case  where  we  obtain  the  minimum  cost,  but  also 
the  maximum  cost  depending  of  the  choice  of  this  angle. 
Similarly,  effects  of  the  Moon’s  eccentricity  and  the  Sun 
are  necessary  to  correct  the  cost  to  maintain  each  rigid 
configuration. 

7.  Dynamical  behavior  of  different  classes  of  solutions 

Using  polar  coordinates,  let  us  define  an  initial 
condition 

V(0)  =  (xA(0)  +  (r(O)cos0(O),r(O)sin0(O)),  vA(0)),  (19) 

as  illustrated  in  Fig.  24.  In  this  last  section,  we  want  to 
study  the  dynamical  behavior  of  the  distance,  denoted  by 
d(Xh(t),X(t)),  between  the  spacecraft  and  the  periodic 


Fig.  24.  Illustration  of  the  relative  position  of  a  spacecraft  around  a 
periodic  orbit  of  Z4. 


orbit  at  each  time  t  using  the  nonlinear  model,  i.e.  without 
linearizing  the  vector  field  A X(t)  =  X(t)  - Xh(t).  First, 
Fig.  25  shows  the  value  of  the  maximum  separation  during 
2  years  for  0°  ^  6(0)  ^  360°  and  r(0)  =  1,2  km  along  the 
long  period  orbit  shown  in  Fig.  2.  Similarly,  Fig.  26  shows 
the  value  of  the  maximum  separation  during  2  years  along 
the  short  period  orbit  shown  in  Fig.  2.  We  have  numeri¬ 
cally  integrated  these  two  initial  conditions  and  the  best 
and  worst  cases  occur  when  6(0)  =  150°  and  6(0)  =  60°, 
respectively.  This  fact  does  not  depend  on  the  initial  condi¬ 
tion  Xh(0).  Therefore,  the  distance  will  be  computed  when 
6(0)  =  150°  and  6(0)  =  60°  to  compare  the  difference 
between  the  quantitative  behavior  of  d(Xh(t),X(t)). 

For  the  first  case,  we  have  taken  four  relative  initial  vec¬ 
tors  r7(0),  r2( 0),  r3( 0),  and  r4( 0)  along  the  direction  of 
6(0)  =  150°  associated  to  the  initial  condition  Xh(0)  of 
the  periodic  orbit  around  L4  (see  Fig.  24).  The  vectors  have 
been  distributed  symmetrically  with  respect  io  Xh(0)  :  r7( 0) 
and  r3( 0)  being  at  an  initial  distance  of  1  km  from  Xh(0) 
and,  r2( 0)  and  r4( 0)  being  at  an  initial  distance  of  2  km. 
We  have  integrated  these  trajectories  during  three  years 
and  the  results  are  shown  in  the  first  row  of  Fig.  27.  As 
we  can  see,  the  maximum  deviation  from  the  starting  sep¬ 
arations  is  less  than  1.5  km  for  the  orbits  with  initial  vec¬ 
tors  rt( 0).  The  distance  d(Xh(t):Xi(t))  is  a  periodic 
function  with  period  equal  to  460  days.  Additionally,  there 
is  no  difference  between  the  qualitative  and  quantitative 
behavior  of  the  distance  for  the  trajectories  starting  at 
the  same  distance  from  2^(0)  :  d(Xh(t),X\(t))  «  d(Xh(t ), 
X3(t))  and  d(Xh(t),X2(t))  «  d(Xh(t),X4(t)). 

Now,  for  the  second  case,  we  take  four  relative  initial 
vectors  q2( 0),  q2(0 ),  q3( 0),  and  q4( 0)  along  the  direction 
of  6(0)  =  60°  corresponding  also  to  the  initial  condition 
2^(0)  of  the  periodic  orbit  around  Z4,  distributed  in  a  sim¬ 
ilar  fashion  to  the  previous  case.  In  the  second  row  of 


0(0) (deg) 


150  155 

e(0)  (deg) 


0(0)  (deg) 

Fig.  25.  Maximum  separation  (in  km)  during  2  years  for  0°  ^  0(0)  <  360°  and  r(0)  =  1,2  km  along  the  long  period  orbit  shown  in  Fig.  2. 
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Fig.  26.  Maximum  separation  (in  km)  during  2  years  for  0°  ^  6(0)  <  360°  and  r(0)  =  1,2  km  along  the  short  period  orbit  shown  in  Fig.  2. 
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Fig.  27.  Distance  d(Xh (t),Xi (t)),  between  the  trajectories  of  r1:  r2l  r5,  r4  (first  row),  and  </7,  q2 ,  q3.  q4  (second  row),  and  the  long  and  short  period  orbits 
shown  in  Fig.  2. 


Fig.  27  the  results  are  shown  for  the  distance 
d{Xh(t),Xi(t)).  Now,  although  the  qualitative  behavior  is 
also  periodic  (period  equal  to  460  days)  and  there  is  no  dif¬ 
ference  between  the  qualitative  and  quantitative  behavior 
of  the  distance  for  the  trajectories  starting  at  the  same  dis¬ 
tance  from  2^(0),  the  deviations  from  the  periodic  orbit  are 
larger  (they  increase  by  a  factor  of  40)  than  the  ones 
obtained  for  the  rt( 0)  initial  vectors  taken  along  the  direc¬ 
tion  of  6(0)  =  150°. 

The  same  study  is  made  along  the  long  and  short  period 
orbits  shown  in  Fig.  7.  The  purpose  is  also  to  analyze  the 
dynamical  behavior  of  the  distance  for  periodic  orbits  that 
are  larger.  In  Fig.  28  it  can  be  seen  that  the  best  position  to 


minimize  the  deviation  is  along  the  direction  of  6(0)  =  149° 
associated  to  the  initial  condition  of  the  long  period  orbit. 
In  the  same  way,  the  position  that  maximize  the  deviation 
is  along  the  direction  of  6(0)  =  59°.  Practically,  the  differ¬ 
ence  with  respect  to  the  previous  results  is  only  1°.  Simi¬ 
larly,  Fig.  29  shows  that  the  position  to  minimize  and 
maximize  the  deviation  is  along  the  direction  of 
6(0)  =  157°  and  6(0)  =  67°,  respectively,  associated  to  the 
initial  condition  of  the  short  period  orbit.  In  this  case, 
the  difference  with  respect  to  the  previous  results  is  7°. 

Similarly,  we  take  four  relative  initial  vectors  along  the 
direction  of  minimum  and  maximum  deviation  corre¬ 
sponding  also  to  the  initial  condition  Xh(0)  of  the  periodic 
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Fig.  28.  Maximum  separation  (in  km)  during  two  years  for  0°  <  6(0)  <  360°  and  r(0)  =  1,2  km  along  the  long  period  orbit  shown  in  Fig.  7. 
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Fig.  29.  Maximum  separation  (in  km)  during  two  years  for  0°  <  0(0)  <  360°  and  r(0)  =  1,2  km  along  the  short  period  orbit  shown  in  Fig.  7. 


orbit  around  Z4,  distributed  in  a  similar  fashion  to  the  pre¬ 
vious  cases.  The  results  are  shown  in  Fig.  30.  In  the  first 
row  we  see  that  the  distance  is  a  periodic  function  (period 
equal  to  460  days  for  the  long  period  orbit  and  545  for  the 
short  period  orbit)  and  that  the  maximum  deviation  from 
the  starting  separations  is  less  than  2  km.  Now,  in  the  sec¬ 
ond  row,  the  distances  are  periodic  (period  equal  to 
460  days  for  the  long  period  orbit  and  545  for  the  short 
period  orbit)  and  the  deviations  from  the  periodic  orbit 
are  also  increased  by  a  factor  of  40.  Additionally,  there  is 
no  difference  between  the  qualitative  and  quantitative 
behavior  of  the  distance  for  the  trajectories  starting  at 
the  same  distance  from  Xh(0). 


Finally,  Figs.  27  and  30  show  that  there  exists  a  linear 
behavior  between  the  distances  for  the  trajectories  starting 
at  1  km  and  2  km  from  Xh(0),  i.e,  the  deviation  for  the  tra¬ 
jectories  starting  at  1  km  is  half  for  the  trajectories  starting 
at  2  km  from  Xh(0).  Remember  that  this  linear  behaviour 
was  also  noted  for  the  cost  to  maintain  a  spacecraft  along 
the  direction  of  zero,  minimum  and  maximum  relative 
radial  acceleration.  Therefore,  there  is  no  difference  in  the 
qualitative  behavior  of  the  distance  for  the  trajectories 
starting  at  1  km  and  2  km  from  Xh(0),  and  the  spacecraft 
that  start  at  those  points  probably  remain  aligned  during 
the  time  span.  To  verify  this  fact,  we  compute  the  difference 
between  the  angles  6 1  and  02  at  each  instant  t  for  all 


F.J.T.  Salazar  et  al.  /  Advances  in  Space  Research  54  (2014)  1838-1857 


1855 


Long  Period  Orbit 

7  f^—r^O)  =  1  km,  a^O)  =  149° 


- r3(0)  =  1  km,  83(0)  =  329" 

- r2(0)  =  2  km,  02<O)=  149° 

---r„(0)  =  2  km,  9^(0)  =  329° 


days 


150 


- q,(0)  =  1  km,  0^0)  =  59° 

—  -q3(0)  =  1  km.  03(O)  =  239° 


0  200  400  600  800 


days 


Short  Period  Orbit 


- ^(0)  =  1  km.  0^0)  =  157° 

- r3(0)  =  1  km.  03<O)  =  337° 

- r2(0)  =  2  km.  e2(0)  =  157“ 


0  200  400  600  800  1000 

days 


140 

120 

100 


- q^O)  =  1  km,  0^0)  =  67° 

___q3(0)  -  1  km.  03(O)  -  247° 

- q2(0)  =  2  km,  02(O)  =  67“ 

- q4(0)  =  2  km,  04(O)  =  247“ 


20 

0 


0  200  400  600  800  1000 

days 


Fig.  30.  Distance  d{Xh{t)1Xi{t)),  between  the  trajectories  of  rl5  r2,  r3,  r4  (first  row),  and  ql5  q2,  q3,  q4  (second  row),  and  the  long  and  short  period  orbits 
shown  in  Fig.  7. 
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Fig.  31.  Difference  between  the  angles  6\  and  02  at  each  instant  t  starting  at  1  km  and  2  km  fromX/j(0)  associated  to  the  initial  conditions  of  the  long  (a) 
and  short  (b)  period  orbits  shown  in  Fig.  2. 


days  days 

(a)  (b) 

Fig.  32.  Difference  between  the  angles  0\  and  02  at  each  instant  t  starting  at  1  km  and  2  km  fromX/j(0)  associated  to  the  initial  conditions  of  the  long  (a) 
and  short  (b)  period  orbits  shown  in  Fig.  7. 
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previous  trajectories  starting  at  1  km  and  2  km  from  A/,  (0). 
In  Fig.  31(a)  and  (b)  we  can  see  the  results  associated  to  the 
long  and  short  period  orbits  shown  in  Fig.  2,  respectively. 
Similarly,  Fig.  32(a)  and  (b)  show  the  results  associated  to 
the  long  and  short  period  orbits  shown  in  Fig.  7,  respec¬ 
tively.  As  we  can  note,  the  difference  of  the  angles  6 1  and 
6 2  at  each  instant  t  is  practically  zero.  Consequently,  the 
spacecraft  remain  aligned  along  the  periodic  orbits  around 
Z4,  but  their  distances  from  the  solution  Xh(t)  oscillate 
periodically. 

8.  Conclusions 

In  this  work  we  have  shown  that,  considering  the  planar 
formation  flight  dynamics  in  the  scenario  of  CRTBP,  Zero 
Relative  Radial  Acceleration  Lines  (ZRRAL)  do  not  exist 
for  periodic  orbits  around  equilateral  equilibrium  points 
when  the  amplitude  of  the  periodic  orbits  is  small  enough. 
For  these  orbits,  the  relative  position  of  the  satellite  has  val¬ 
ues  for  the  angle  6  such  that  its  relative  acceleration  has 
either  maximum  or  minimum  radial  component.  On  the 
other  hand,  the  ZRRAL  exist  when  the  amplitude  of  the 
periodic  orbits  is  large  enough.  However,  for  these  orbits, 
the  ZRRAL  do  not  exist  for  every  point  of  the  periodic 
orbit.  Therefore,  in  these  cases,  the  relative  position  of  the 
satellite  has  values  for  the  angle  0  such  that  its  relative  accel¬ 
eration  has  maximum,  minimum  or  zero  radial  component. 

The  cost  estimate  to  maintain  a  satellite  following  an 
artificial  trajectory  has  been  determined  by  the  integral  of 
the  residual  acceleration.  In  the  case  of  periodic  orbits 
around  Z4,  it  has  been  shown  that  the  cost  estimate  is  less 
when  the  relative  position  of  the  satellite  points  towards 
the  direction  where  the  radial  component  of  the  relative 
acceleration  is  minimum,  although  the  ZRRAL  do  exist. 
Due  to  the  fact  that  these  regions  appear  and  disappear 
along  the  periodic  orbit,  producing  a  discontinuity  in  the 
components  of  the  relative  velocity  and  acceleration  and, 
therefore,  a  higher  cost.  On  the  other  hand,  the  cost  does 
not  depend  on  the  amplitude  of  the  periodic  orbit  when 
the  spacecraft  follows  the  direction  where  the  relative 
radial  acceleration  component  is  maximum.  However,  the 
cost  increases  with  the  amplitude  when  the  spacecraft  fol¬ 
lows  the  direction  where  the  radial  acceleration  component 
is  minimum  or  zero.  Similarly,  we  have  computed  the  cost 
estimate  to  maintain  a  satellite  when  the  relative  position 
was  fixed,  i.e.  when  the  angle  6  remains  constant  along  a 
family  of  periodic  orbits  around  Z4.  The  function  cost 
has  a  maximum  at  0  &  60°  and  a  minimum  at  0  «  150°. 

A  formation  of  three  satellites  has  been  studied.  The 
configuration  is  a  triangular  geometry,  such  that  the  rela¬ 
tive  position  along  the  periodic  orbit  is  fixed  and  forms 
an  equilateral  triangle.  The  cost  estimate  to  maintain  this 
configuration  is  minimum  when  one  of  the  satellites 
remains  on  the  periodic  orbit. 

Note  that  in  the  force  model  defined  by  the  CRTBP,  the 
cost  estimate  to  maintain  the  satellites  in  these  configura¬ 
tions  (ZRRAL,  maximum  or  minimum  radial  acceleration, 


fixed  angle)  is  very  low.  This  fact  is  due  to  the  stability  of 
the  periodic  orbits  around  Z4,  so  that,  particles  that  begin 
very  close  to  a  periodic  orbit,  their  mutual  distances  will 
keep  practically  fixed,  so  the  control  applied  to  the  space¬ 
craft  is  very  small.  However,  perturbations  like  solar  radi¬ 
ation  pressure,  solar  gravitational  attraction  and  Moon’s 
eccentricity  affect  the  stability  of  the  periodic  orbits  around 
Z4.  Therefore,  in  subsequent  works,  all  this  kind  of  pertur¬ 
bations  will  be  considered  to  compute  a  more  accurate 
cost. 

The  dynamical  behavior  of  the  distance  between  a  con¬ 
stellation  of  satellites  and  the  periodic  orbits  around  L4  has 
also  been  determined.  It  can  be  seen  that  the  maximum  sep¬ 
aration  function  does  not  depend  on  the  initial  condition 
taken  along  the  periodic  orbits  and  that  this  function  has 
a  period  of  180°.  Similarly,  the  best  position  to  reduce 
the  maximum  separation  is  about  6  =  150°.  On  the  other 
hand,  the  position  that  increase  the  maximum  separation 
is  about  6  =  60°.  Thus,  placing  satellites  in  the  directions 
of  either  minimum  or  maximum  deviation  from  the  start¬ 
ing  separations,  the  distance  from  the  periodic  orbit 
behaves  as  a  periodic  function  in  such  a  way  that  the  satel¬ 
lites  remain  aligned.  Therefore,  if  we  want  to  maintain 
fixed  the  configuration,  a  thrust  would  be  necessary  only 
to  control  the  separation.  In  the  case  of  minimum  devia¬ 
tion,  since  the  deviation  in  some  cases  was  less  than 
2  km,  the  cost  estimate  to  maintain  the  distance  fixed 
would  be  very  small.  Therefore,  simulations  have  shown 
that  the  direction  of  minimum  radial  component  defines 
the  initial  conditions  for  the  satellites  with  respect  to  the 
periodic  orbit  around  Z4  so  that  the  drift  between  them 
be  minimized.  This  information  has  a  practical  value  for 
the  control  of  the  formation  flying  since  the  fuel  spent  will 
be  less  if  the  spacecraft  are  set  up  in  the  direction  where  the 
relative  radial  acceleration  component  is  minimum. 
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